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J> I Abstract 
OO 

i We discuss all contributions from the Zweig-rule-satisfying SU(3)-breaking 

On I final state interactions (FSIs) in the B PP decays (neglecting charmed 

I intermediate states), where PP=Tnr, ttK, KK, ■Kri(T]'), and Kr]{r]'). First, 

I effects of SU(3) breaking in rescattering through Pomeron exchange are stud- 

i-C I ied. Then, after making a plausible assumption concerning the pattern of 

I I SU(3) breaking in non-Pomeron FSIs, we give general formulas for how the 

latter modify short-distance (SD) amplitudes. In the SU(3) limit, these for- 
' mulas depend on three effective parameters characterizing the strength of all 

non-Pomeron rescattering effects. We point out that the experimental bounds 
^ ■ on the B — > K'^K^ branching ratio may limit the value of only one of these 

' FSI parameters. Thus, the smallness of the B K^K^ decay rate does not 

imply negligible rescattering effects in other decays. Assuming a vanishing 
value of this parameter, we perform various fits to the available B — > PP 
branching ratios. The fits determine the quark-diagram SD amplitudes, the 
two remaining FSI parameters and the weak angle 7. While the set of all 
B — > PP branching ratios is well described with 7 around its expected Stan- 
dard Model (SM) value, the fits permit other values of 7 as well. For a couple 
of such good fits, we predict asymmetries for the B — > Xvr, TT^ri{r]'), K^rj{'q') 
decays as well as the values of the CP-violating parameters Stttt and Ct^t^ for 
the time-dependent rate of B^{t) vr+vr". Apart from a problem with the 
recent -B"*" TT~^r] asymmetry measurement, comparison with the data seems 
to favour the values of 7 in accordance with SM expectations. 

PACS numbers: 13.25.Hw, 11.30.Hv,12.15.Hh,11.80.Gw 
* E-mail: zenczyko@iblis.ifj.edu.pl 
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1 Introduction 



The majority of the analyses of CP-violating effects in B decays assume that the 
relevant amplitudes are given by short-distance (SD) expressions only. In particular, 
for B decays into two pseudoscalar mesons {B —>■ PP), any possible final state 
interactions (FSIs) are usually completely neglected. It is very difficult to assess if 
this neglect is justified or not. Some authors have argued that such effects should be 
negligible [1, 2] since the B mass is already quite large. In other papers it is stressed 
that the FSIs should be important and that any reliable analyses of B decays must 
take these interactions into account [3, 4, 5, 6, 7]. It has been suspected that the 
inelastic FSIs are particularly important [3, 6]. Unfortunately, with our insufficient 
knowledge of the PP interactions at 5.2 GeV, there is virtually no hope that the 
relevant rescattering effects may be calculated reliably. 

In order to overcome this obstacle, in a recent paper [8] we analysed an SU(3)- 
symmetric approach with the built-in Zweig rule, in which our ignorance as to the 
size of inelastic rescattering was reduced to a set of only three effective (complex) 
parameters jointly describing all inelastic final state interaction (IFSI) effects. It 
was shown that the SU(3)-symmetric rescattering leads to a simple redefinition 
of quark-diagram amplitudes, thus permitting the use of a diagram description in 
which, however, weak phases may enter in a modified way. Furthermore, a simple 
estimate was made as to the size of error which could be committed while extracting 
the value of the unitarity-triangle angle 7 when such modifications are not taken 
into account. 

In the present paper, we extend the general scheme of ref.[8] and introduce 
SU(3) breaking both in the elastic and in the inelastic final state interactions. The 
introduction of SU(3) breaking makes it reasonable to attempt a detailed description 
of the data. When doing so, we take into account all short-distance amplitudes 
usually considered as the dominant ones (Section 2), and make certain assumptions 
as to the form of FSIs and SU(3) breaking (Sections 3 and 4). In Section 4 we 
also discuss at some length the point that estimating the size of all rescattering 
effects on the basis of the B — > KK data is significantly more difficult than usually 
acknowledged. Then, in Section 5, we perform fits to the experimental branching 
ratios of the B — PP decays, and discuss their implications. A brief summary 
appears in Section 6. 
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2 Short-distance amplitudes 



Short-distance amplitudes may lead not only to the PP states but also to the general 
M1M2 states, with Mj representing various heavy mesons. Consequently, the PP 
pair observed in B decay may be produced in three ways: it may not participate 
in any rescattering after being produced in a SD process, it may undergo elastic 
rescattering, and, finally, it may result from inelastic rescattering of M1M2 into PP. 
As discussed in [8], with the help of the unitarity condition, contributions from other 
inelastic intermediate states (such as many-body states M\M2...Mn) may be always 
incorporated into the contribution from M1M2. 

All SD amphtudes B M1M2 may be classified in the same way as stan- 
dard SD amplitudes B PP, ie. T.T' (tree), C.C (colour suppressed), P,P' 
(penguin), E,E' (exchange). A, A' (annihilation), PA, PA' (penguin annihilation), 
S,S' (singlet penguin), SS,SS' (double singlet penguin). As usual, we denote 
strangeness-conserving AS* = (strangeness-violating \AS\ = 1) processes by un- 
primed (primed) amplitudes respectively. Electroweak penguin contributions may 
be included via the replacements: T — ^ T + P^y^/, P P — P|^y/3, C ^ C + Pew, 
S ^ S — Pew/'^ [9] (with analogous expressions for the primed amplitudes). 

The essential assumption of refs.[10, 8] is that the tree, penguin, etc. amplitudes 
for the production of various M1M2 states are proportional to the corresponding 
amplitudes for the production of the PP pair. One may argue that the relevant 
coefficients of proportionality are approximately independent of the diagram type 
(tree, penguin, etc.) considered. The common remaining single coefficient of propor- 
tionality may be absorbed into the rescattering amplitudes M1M2 PP, for which 
the Zweig-rule is assumed. Finally, the sum over all intermediate states M1M2 may 
be performed leading to the appearance of only three effective complex parameters 
representing the relevant sums and corresponding to the presence of three Zweig-rule 
satisfying SU(3)-symmetric forms for M1M2 — > PP (for more details, see [10, 8]). 

As a result of these simplifications, all contributions from various short-distance 
B M1M2 amplitudes get expressed in terms of relevant standard B — > PP short- 
distance amplitudes. Our whole approach to inelastic rescattering depends therefore 
on standard T, P, ...P', T' , ... etc. amplitudes (with appropriate weak phases) and on 
parameters effectively describing the rescattering. In order to simplify the discussion 
and study the effect of FSIs only, we assume that the strong SD phases are negligible. 
(In ref.[l] these phases were estimated to be of the order of 10°, while in ref. [11] 
it is argued that the FSI-uncorrected "bare" amplitudes do not contain any strong 
phases - see the comment after Eq.(16) therein). This assumption may be relaxed 
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in future. 



Some of the SD quark-diagram amplitudes are related. In an approach in which 
FSIs break SU(3), one should incorporate SU(3) breaking into the SD relationships 
as well. Therefore, we assume that the tree SD amphtudes satisfy the following 
relation [12]: 

V f 

T' = ——Tsa 0.276 T (1) 

Vud /tt 

Both tree amphtudes have the same (weak) phase: T/\T\ — T'/\T'\ — e*'"'. 

The penguin SD amplitudes are dominated by the t quark, so that the weak 
phase factor is e"*^ for P and ±1 for P' (ie. P' — ±\P'\). We use the estimate [12] 



P = -e-^^ 



td 



ts 



P' ^ -0.176 e-'^P'. (2) 



In the fits of Section 5, we accept (3 = 24°, which is in agreement with the world 
average [13] sin 2/3 = 0.734 ±0.054. We accept (as it is usually done) that the value 
of the penguin SD amplitudes docs not depend on the flavour of the quark- ant iquark 
pair created to produce the M1M2 state. For example, standard SD contributions 
from penguin P in i?^ — > tt+tt^ (or tt^tt"), and in 5+ K^K^ are given by SU(3) 
considerations only, despite the fact that in these two processes the produced quark- 
antiquark pairs are of different flavours. 

We accept the relations between the tree and the colour-suppressed amplitudes 
given by the SD estimates: 

C^iT (3) 

and 

c' = r'(e-(i + OW-^^) (4) 

where we take ^ = 0.17, assuming C, ~ 0.42, ie. midway between 1/Nc and 

the value of 0.5 suggested by experiment, and using Ci ^ —0.31 and C2 ~ 1.14 [14]. 
The contribution from the electroweak penguin P^w been included in Eq.(4), 
with 5ew ~ +0.65 [15] (other electroweak penguins are neglected). 

The last independent SD amphtude considered here is the singlet penguin amph- 
tude 5", whose weak phase is (data requires that this amplitude be sizable [16, 12]). 
Thus, the SD amplitudes and our whole approach depend on four SD parameters: 
\T\, P\ 5", and the weak phase 7. The remaining SD amplitudes {E, E' , S, PA, ...) 
are assumed to be negligible. 
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3 SU(3)-breaking in Pomeron-exchange-induced 
rescattering 



If we gather all SD amplitudes B PP (as well as those of B ^ M1M2) into vector 
w, and accept that FSIs cannot modify the probabilities of the original SD weak 
decays, it follows that vector W representing the set of all FSI-corrected amplitudes 
is related to w through [6, 10]: 

W = S^/V, (5) 
(in the one-channel case, Eq.(5) reduces to the Watson's theorem [17]). 

Let us consider now elastic PP rescattering only (ie. with w restricted to its 
part corresponding to B ^ PP processes, and similarly for W). For high energies 
this rescattering is approximately independent of energy. We shall use Regge ter- 
minology and call this energy-independent term a Pomeron-induced contribution. 
Since Pomeron exchange is known to be substantial, the B PP amplitudes at 
s = m\ should be corrected for Pomeron-induced rescattering. Treating Pomeron- 
induced FSIs as a small correction to the SD expressions for B —>■ PP amplitudes 
corresponds to expanding S^/^ = (1 + iT)^/^ = (1 + 2iAy/^ = 1 + iA + ... and 
keeping terms linear in A only. Thus, one gets [6]: 

W^{l + iA)w (6) 

Because the amplitudes for Pomeron exchange are predominantly imaginary, we 
have 

A = ia (7) 

with real a. In the SU(3)-symmetric world, all elements of a are identical, and their 
common value is a ~ 0.16 (cf [3] and Eqs.(10,17) in [6]). Consequently, Pomeron- 
induced rescattering rescales all SD amplitudes in the same way: 

W = (1 - a)w. (8) 

It is only when SU(3) is broken that the rescaling is different for different decay 
channels, and deviations from the standard SD form could be observed in principle. 

When SU(3) is broken, the values of a differ for different final channels P1P2. In 
a simple model for Pomeron used in [3, 18], they are given by 



5 



with the values of /?7r, Pk (meson- Pomeron couphngs) and ^tt, (slope coefficients 
for the relevant couplings) extracted from data on Trp and Kp scattering. In the 
following we will use the averages of values given in [3, 18], ie.: 

(5^ = 3.47 A/mb 

Pk = 2.78 

K = 1-93 GeV"^ 

hK = 0.9 GeV-2 (10) 

In order to estimate and 6^, 6^/, we assume perfect mixing for 77,77' (ie. 

77 = {uu + dd — ss)/y/3, and 77' = {uu + dd + 2ss)/y/6) corresponding to the octet- 
singlet mixing angle of ^ = —19.5°, (see eg. [19, 20, 21, 16]; for a different approach 
to 77 — 77' mixing in S — > ^77*^') decays see [22]), and derive [18]: 

= (p^ + 2/3j^)/3 ^ 3.01 

Pv' = i-P^ + 4/9i^)/3 ~ 2.55 Vmb 

= {b^ + 2bK)/^ ^ 1.24 

b^, = (-6^ + 46^)/3fa 0.56 GeV-2 (11) 

Note that for the Kfj' channel the denominator in Eq.(9) is particularly small. In 
this channel the Pomeron-exchange-induced correction is therefore relatively large 
which may possibly affect the extraction of the short-distance S' amplitude from 
the data. 

The resulting pattern of SD amplitudes corrected for Pomeron-induccd rcscat- 
tering differs from standard SD expressions by departures from SU(3) only. Con- 
sequently, we introduce SU(3)-symmetric rescaled amplitudes T, T', P, P', ..etc., 
defined as 

fC) = r(')(l-a(7r7r)) 

pC) = P(')(l-a(7r7r)) (12) 

, 

and the SU(3)-breaking corrections K{PiP2) = {a{7in) — a(PiP2))/(l — a(7r7r)). The 
complete set of SD amplitudes corrected for SU(3)-breaking Pomeron-exchange- 
induced rescattering is gathered in Table 1. 
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Table 1: SD amplitudes for decays B^,B^ — > P1P2 corrected for SU(3)-breaking 
Pomeron-exchange- induced rescattering 



Decay 


rescaled and corrected SD 










t\ t\ 


—r'yl + KyKK)) 






1 CT J- r' J- J- Ki'w'nW 

T^U + + '^^A^ + -"^ 1"^^?;; 








dO 
Bd 


— » TT^TT 


-{f + P) 



















-P{1 + K{KK)) 


B+ 




-P'{1 + K{TrK)) 






-^{f + C' + P'){l + K{7rK)) 




r]K+ 


-j^_if + C' + S')il + KirjK)) 




r]'K+ 


^{T + €' + 3P' + 4S')(1 + K{iK)) 




n-K+ 


{f + P'){1 + K{ttK)) 






j.{C'-P'){l + K{7rK)) 












-^{C + 3P' + 45') (1 + K{r]'K)) 



4 Inelastic SU (3) -breaking FSI with Zweig rule 

Analysis of inelastic SU(3)-breaking effects follows the approach of [8]. As in ref. [8], 
in the present paper we do not consider contributions from intermediate charmed 
states (thus neglecting the long-distance "charming penguins"). Since they may 
be important [23, 24, 25, 26, 11], their analysis merits further work. The most 
general Zweig- rule-satisfying rescattering M1M2 P1P2 is described by two types 
of connected diagrams: the "uncrossed" diagrams of Fig.l(M), and the "crossed" 
diagrams of Fig.l(c). By virtue of Bose statistics, the final P1P2 pair must be 
in an overall symmetric state. Our definition of inelastic rescattering includes a 
non-Pomeron contribution from P1P2 — > P1P2 transitions, which - together with 
the Pomeron-exchange-induced part of these transitions - are usually classified as 
elastic. 
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4.1 SU(3)-invariant rescattering amplitudes 



In the SU(3) case, the requirement of Bose statistics for P1P2 means that there are 
two types of uncrossed M1M2 -P1-P2 amphtudes, ie. (using a particle symbol for 
the corresponding SU(3) matrix): 

Tr{{MlMl}{P^,P2})u+ (13) 

and 

Tr([M|,Ml]{Pi,P2})H- (14) 

where the requirement in question is reflected by the presence of the anticommutator 
{Pi, P2} of meson matrices, and u± denote the strength of rescattering amplitudes. 
Eqs. (13,14) incorporate nonet symmetry for both intermediate and final mesons. As 
explained in [8], invariance of strong interactions under charge conjugation demands 
that mesons Mi and M2 belong to multiplets of the same (opposite) C-parities for 
the first (second) amplitude above. 

For the crossed diagrams, the requirement of Pi # P2 symmetry admits one 
combination only [8]: 

Tr(Mi^PiM|P2 + Mi^PaMjPi) c (15) 

where c denotes amplitude strength. This combination is symmetric under Mi # M2 
as well. Consequently, it is charge-conjugation invariant if Mi and M2 have C- 
parities of the same sign. 

4.2 Modifications due to SU(3) breaking 

We will incorporate SU(3) breaking into the FSl amplitudes of Eqs. (13, 14, 15) in 
the simplest possible way. First let us consider -u-type diagrams (Fig.l(u)). In these 
diagrams one quark (or antiquark) from meson Mi ends up in the final pseudoscalar 
meson, while the other one annihilates an antiquark (quark) from meson M2. It is 
well known that such quark-antiquark annihilations are suppressed when the rele- 
vant qq pair has high energy, and that they are suppressed even more strongly for 
the ss pair. In the Regge language, the first statement corresponds to meson ex- 
changes being suppressed at high energies, the latter - to the fact that intercepts of 
Regge trajectories for mesons containing strange quarks lie below those for mesons 
composed of u, d, u, d only. The additional suppression of ss annihilation with re- 
spect to that of uu (or dd) depends on the energy of the qq pair. Since we want 
to analyse the main effect of SU(3) breaking only, we assume that an exchange of 
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a strange (anti)quark between mesons Mi and M2 (or between Pi and P2) is sup- 
pressed by the same factor (e) for all intermediate states. On the other hand, the 
amplitudes for the uncrossed diagrams in which strange (anti) quarks from mesons 
Ml end up in final pseudoscalar mesons (ie. are not annihilated) are not suppressed 
by SU(3)-breaking effects. 

The relevant u-type amplitudes may be then calculated from the appropriate 
generahzations of Eqs. (13,14). For the contribution from mesons Mi and M2 of the 
same charge-conjugation parities (C(Mi)C(M2) = -|-1) we have, for example: 

ITtHmUM + MlhMl){PihP2 + P2hPi)) u+ 
+ iTr((Mr/,Mf + Ml^mDiPllePl + PlleP^) u+ (16) 



where 



1 
1 
e 



(17) 



In Eq.(16) we divided the whole contribution into two parts, depending on whether it 
is the strange quark or antiquark from (say) Mi which is annihilated. Contributions 
from the C(A'/i)C(M2) = —1 states may be calculated in a similar way (one has to 
remember that the negative sign between MiI^M2 and M2leMi is cancelled by the 
negative sign in the (antisymmetric) wave function of C(Mi)C(M2) = —1 two-meson 
states) . 

Since SU(3) is to be broken, the choice of definite SU(3) representations for the 
intermediate M1M2 states is not appropriate. Admitting the linear combinations of 
27, 8s, 8{8i}, 1{88}) ^'^d 8a (considered in [8]) is not sufficient cither, since 

for broken SU(3) the complete set of M1M2 intermediate states contains admixtures 
from other SU(3) representations. If all the C(Mi)C(M2) = ±1 intermediate states 
are to be taken into account properly, one may first list all states of definite charge, 
strangeness and isospin, and composed of two mesons of definite type, ie. with 
fiavour quantum numbers of nK orrjK or... . These states may be ordered (in the 
sense that ttK is different from Kir) or, alternatively, their symmetric or antisym- 
metric combinations (under n K etc. interchanges) may be formed. Then, SD 
decay amplitudes into these states have to be evaluated. Finally, the sum over the 
contributions from all such states has to be carried out. 

We have performed all the necessary calculations with the result that the sum 
over all C{Mi)C{M2) = ±1 intermediate states leads to the formulas given in the 
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second column of Table 2, where 



1 



u+ + 1 



u 



u 



1 — a(7r7r) 
1 



2 1 - a(7r7r) 



d 



d 



1 — a(7r7r) 



(m+ — - 

1 — a(7r7rj 



(18) 



and 



A' = 



A 



((2 + e)P + f) d 

{{2 + e)P' + f' + eS') d. 



(19) 



Thus, the inelastic corrections are given in terms of the products of the SD ampli- 
tudes and the FSI parameters (here: u, d). For example, there may be a contribution 
proportional to Td. Since we finally express all formulas in terms of the amplitudes 
modified for Pomeron-induced rescattering (eg. in terms of T = T(l — a(7r7r)) etc.), 
in Eq.(18) we introduced the rescaled FSI parameters u and d so that eg. Td — Td. 
For completeness, in Table 2 we give formulas for the S° decays as well. 

We incorporate SU(3) breaking into the c-type amphtudes in a completely anal- 
ogous fashion. Namely, we assume that strange (anti) quark interchanges are sup- 
pressed by factor e (in general, this factor may be different from that used for it-type 
diagrams) . The relevant c-type amplitudes may be then calculated from an appropri- 
ate generalization of Eq.(15). As pointed out in [8], charge conjugation invariance 
of strong interactions requires that only symmetric M1M2 states contribute. For 
broken SU(3), Eq.(15) is replaced by 



As in Eq.(16), we divided the whole contribution into two parts depending on 
whether it is the strange quark or antiquark from (say) Mi which is exchanged. 
Using the above expression and the expressions for the SD amplitudes, and sum- 
ming over all the intermediate states, one obtains the corrections induced by the 
c-type IFSIs. They are listed in the third column of Table 2, where 



In the hmit of e ^ 1, all formulas of Table 2 reduce to those given in [8], while for 
e = SU(3) is maximally broken. 



+ \l^{MrhP^MthPl-rMri,P^MrhPl) c 



(20) 



c = c/(l — a(7r7r)) 



(21) 
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Table 2: Inelastic SU(3)-breaking rescattering contributions: A = ((2 + e)P + T) d; 
A' = ((2 + e)P' + f + eS') d 



Decay uncrossed FSI diagrams crossed FSI diagrams 



B+ - 







V 




K+K^ 


-e(A + 2uC) 









__|(A + 2zZC) 


-^2c{T + C + P{2-e)) 






_^(A + 2SC) 


_ 1 2c(r+C' + 2P(l + e)) 






-(A + 2u(T + 2P)) 


-2cC 






^(A + 2n(f + 2P)) 






K+K^ 


2u{eT + (1 + e)P) 







K°K° 


-eA-2w(l + e)P 









-^(A-2zZT) 








-^{A-2uf) 






^ n+K- 


-A 


-(l + e)c(7 








-^(l + e)c-f 








-^c((2-e)P + T) 




r]'K° 




-^c(2(l + e)P + r) 


B+ - 


n+K° 


-A' - 2u(C" + S') 


-c{l + e)S' 






^(A' + 2«(C' + y)) 


^c(l + e)(r + C' + 5') 




r]K+ 


^(l-e)(A' + 2n(C" + 5')) 


^c(l + e)(r' + C' + P'(2 - e) + ^(l - e) 




r]'K+ 


i±|^(A' + 2«((7' + 5')) 


^c(l + e)(T' + C' + 2(1 + e)P' + (1 + 2e)5') 


Bl- 


^ ■K-K+ 


A' + 2uS' 


c(l + e)(C' + 5') 






-^(A' + 2S^') 


^c(l + .)(f'-y) 






^(l-e)(A' + 2^5') 


^c(l + 6)(f + (2-e)P' + (l-6)5') 




r]'K° 


l^iA' + 2uS') 


^c(l + e)(r + 2(l + e)P' + (l + 2e)y) 


Bl- 


-> 7r+7r~ 


-2eu{2P' + T') 
V2eu{2P' + f') 








K+K- 


A' + 2eu((l + e)P' + eT' + 5') 


2ec(C" + 5') 




K°K° 


-A' - 2eu((l + e)P' + S') 


-2ec5' 




7r% 
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4.3 Structure of full FSIs 



For small inelastic final-state interactions, Eqs.(5,6) suggest the following approxi- 
mation of all FSI effects: 

W w - aw + iAWinei (22) 

where the three terms correspond to the contributions from the unmodified SD 
amphtudes, the Pomeron-exchange-induced corrections, and the inelastic FSI cor- 
rections (including the P1P2 — > P1P2 elastic transitions not mediated by Pomeron) 
respectively. Here AWinei (proportional to J2miM2'^\^i^2){MiM2\'w ) is given 
by expressions for the inelastic FSIs gathered in Table 2. For negligible strong SD 
phases, it is the third term in Eq.(22) which allows the existence of direct CP vio- 
lation effects. This term provides a specific prescription for how strong phases are 
generated by quark interchanges between outgoing mesons. In other words, the pat- 
tern of FSI phases in all B — ^ PP decays is governed by three (in general complex) 
parameters d, u, c corresponding to different fiavour-fiow rescattering topologies and 
by the value of the SU(3)-breaking parameter (s) e. 

4.4 Size of rescattering effects and B — ^ KK decays 

The Pomeron-induced FSIs and a contribution from non-Pomeron-mediated transi- 
tions P1P2 P1P2 together comprise elastic rescattering. The non-Pomeron contri- 
butions to elastic transitions (eg. quark-line exchange diagrams for tt+tt" tt+tt^) 
should be treated alongside symmetry-related contributions (ie. tt+tt^ 7i"7i" or 
Tr^n" K^K^ etc.), as they all have common origin. For the SU(3) case all such 
"quasi-elastic" P[P2 P1P2 transitions were estimated in the Regge approach [27]. 
The resulting differences between strong phases in the singlet, octet, and 27-plet PP 
channels (see also [18]) vanish at high energy, while at the S-meson mass they turn 
out to be nonnegligible yet small, of the order of 10°. Consequently, inclusion of full 
elastic FSIs should not lead to a significant change in the quality of data description 
(see also the fits of the next Section) . 

As for the inelastic rescattering. Table 2 provides the basis for the relevant dis- 
cussion. 

If FSIs satisfy SU(3) (ie. if e = 1), all the A and A terms in Table 2 may be 
absorbed into the new redefined penguins [8] (compare Eq.(22)): 

P = P + iA 

P' = P' + iA (23) 
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With our assumptions of SU(3)-symmetric SD penguin amplitudes (cf. comment 
after Eq.(2)), such a redefinition is possible only if e = 1 (compare the relevant 
A-dependent corrections to K^K^ and ti^tt~ in Table 2). As can 

be seen from the presence of the SD tree amplitudes in the redefined penguins in 
Eq.(23) (cf. A = (3P + T)(i), parameter d is related to the size of the long-distance 
(■u-quark-loop) penguin. Formulas (23) indicate that contributions from penguin 
topologies with internal w-quark loops cannot lead to significant modification of 
total amplitudes - all such effects are consistently absorbed everywhere into new 
redefined penguins P, P'. The only change is in the phase factors since A's include 
terms depending on 7. In general this leads to nonzero asymmetries, and should 
affect the determination of 7, as the (effective) amplitudes will now interfere in a 
different way. 

In some papers it was argued that the B K^K~ decays could provide an 
estimate of the size of rescattering effects. Note, however, that this decay amplitude 
depends on u only. The B — > K'^K~ branching ratio is independent of d and 
c, and, consequently, the size of long-distance penguins is not restricted by 5 — > 
K'^K" . This means that this branching ratio is not such a good place to estimate 
the "typical" size of FSI effects as it has been thought so far. 

It is also sometimes said that the size of rescattering effects may be gleamed from 
the B^ K~^K^ decay which is related to the B^ — » ir'^K^ decay by an interchange 
of all down and strange quarks [29]. Here the standard argument assumes U-spin 
flavour symmetry of strong interactions. When SU(3) is broken, a look at Table 2 
and Eqs. (19,23) shows that the conclusions from the comparison of S+ — > K^K^ 
and B^ — > tt+X^ decays cannot be obtained in such a simple way as originally 
thought. Namely, with the contribution from fZ-generated FSI effects being bounded 
by the smallness of the B^ — > K^K' branching ratios, the FSI effects in S+ — > 
K'^K^ are proportional to term eA. However, on the basis of Regge ideas and 
our knowledge of high-energy multiparticle production processes in which KK pairs 
are rarely produced, one expects that e is small. (The assumption of negligible e 
seems to be corroborated by the e-dependence of our fits below.) Consequently, the 
rescattering term in B~^ K~^K^ could be smaller by a factor of e from what is 
expected on the basis of U-spin symmetry with B^ — ^ n^K^. Therefore, despite 
the relative 1/A^ factor [29], the overall FSI effects in B'^ — > K~^K^ need not be 
much larger than those in B^ n^K^. Thus, from the smallness of FSI effects in 
B^ — i> K^K^ one cannot infer that such effects are negligible elsewhere. In fact, a 
A-induced term, such as that in S+ — > K'^K^, is present in all formulas in which 
the SD penguin P contributes. This leads (in the SU(3) hmit) to the replacement 
of the original SD penguin amplitude P by the effective penguin P given by Eq. 
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(23). It is only through a combined description of all the B PP branching 
ratios (and possibly asymmetries) such as these attempted in this paper (ic. not 
just of — i> K~^K^ and B^ — > n^K^ decays) that the effects induced by terms 
proportional to A can be hopefully determined. 

In order to study only the most important effects, we make now three assump- 
tions for the fits of the next Section: 

1) First, we put e = thus breaking SU(3) maximally. 

2) Second, the present upper bound on the value of the B — > K^K~ branching 
ratio (< 0.6 • 10~^) limits the size of u quite severely. Thus, we assume for simplicity 
that M = 0. 

3) Third, with no bounds set hy B ^ K^K^ on d and c we must treat these 
parameters as free. However, while the value of d could be complex, one expects 
that c should be real (as required by the condition of no exotics in the s-channel - 
see Fig. 1(c); for the Regge model the corresponding expressions may be found in 
ref.[27]). Consequently, we will have three real parameters: Re (J, Im d. Re c. 

With e fixed, our formulas depend on six real parameters: |T|, P', S' , Red, Imd, 
Rec (in addition to weak phases). This may be compared to the approach of [28] 
which is less specific as to the origin of strong phases and involves seven independent 
hadronic parameters. 

5 Fits 

In order to estimate the effects which SU(3)-breaking rcscattcring may induce, we 
performed fits to the available branching ratios of B decays. We decided to compare 
the case with no FSIs (or with SU(3)-symmetric Pomeron-induced FSIs only) to the 
following two cases: 

(a) SU(3)-brcaking Pomeron-exchange-induced FSIs only, 

(b) both Pomeron-exchange-induced and non-Pomeron SU(3)-breaking FSIs. 
Since one of the objectives of this paper was to test the FSI effects, we assumed 
that the relative strong SD phases are negligible, ie. all direct CP violation effects 
involve only the long-distance strong phases generated by the FSI term iAWi„ei in 
Eq.(22). 

As the data constraining our fits we used only the branching ratios of the B — >■ 
PP decays (ie. we did not inchide the data on asymmetries). The first and the 
second column of Table 3 specify the decay channels i considered and the values 
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of the experimental branching ratios (and their errors) as used in our paper (from 
[12, 30, 31]). In the calculations themselves, the branching ratios were corrected for 
the deviation of the ratio of the tb+ and tbq lifetimes from unity. The sum over all 
these decay channels i of the deviations between the experimental and theoretical 
branching ratios Bj normalized to their experimental errors: 

nSD amp,; FSI ^™™.) = ^ . (24) 

was subject to the minimization procedure (see eg.[l, 32]). Note that in our fits 
we used not only the S — > tttt and B — > ttK branching ratios (as in [1]), but also 
the remaining B PP branching ratios not considered elsewhere (in particular, 
those for B — > Kr],Kr]'). Wc performed several different fits, first keeping some of 
the arguments of / in Eq.(24) fixed, and then letting them free. The minimization 
procedure gave the best values of T, P', and S' (for different values of weak phase 
7) as well as the values of the FSI parameters. The fits permitted predictions of 
CP asymmetries in i? — > Kn, the values of parameters S'tttt and C,r7r describing the 
behaviour of the time-dependent rates in B2{t) — > 7t~^7t~ , etc. Below we discuss our 
results in more detail. 



5.1 Pomeron- induced rescattering 

Consider first the situation with Pomeron- induced FSIs only, ie. d = c = u = 0. 
Two cases differing with respect to the sign of (real) P' may be distinguished. A 
negative value of P' corresponds to vanishing differences of SD strong phases (eg. 
5pi, St'), while its positive value corresponds to this difference being 180°. Using P', 
\T\, and S' as free parameters, we minimized / for different values of 7 for the no-FSI 
case (all a's vanish), and for case (a) above. Dependence of the minimum value of / 
on the value of 7 is shown in Fig. 2. From Fig.2 (and Table 3) one can see that the 
introduction of SU(3)-breaking Pomeron-induced FSIs does not lead to a significant 
improvement in the description of data. Since non-Pomeron contributions to elastic 
rescattering cannot be large at B mass, this result is in contradiction with a recent 
paper [33] which claimed that data provide evidence for a large effect due to SU(3) 
breaking in elastic rescattering. 

The preferred values of 7 are in the range of around 85^ < 7 < 125° (0° < 7 < 
60°) for P' < (P' > 0). The best fit is obtained for P' < with 7 ^ 102° (see Table 
3), in agreement with earher determinations preferring 7 > 90° [1, 34, 35]. Such a 
large value of 7 is in disagreement with the estimates in the Standard Model, which 
lead to 75M ~ 64.5° ±7° [34], or, more conservatively, to the region of 50° < 7 < 80° 
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(see eg. [36, 37, 38]). The approach of ref.[28] permits shghtly smaller values of 7, in 
the range of approximately 75° — 85", at the cost of introducing seven independent 
parameters in place of |T|, P', and S' (see also the next subsection). Table 3 shows 
that the inclusion of SU(3)-breaking Pomeron-induced FSIs enhances the value of 
the S'/P' ratio when extracting it from data. 

5.2 Inelastic rescattering 

Since even when u — e — there are still three real FSI parameters (Re J, Im J, 
Re c), it is instructive to consider first the two limiting cases when 1) \d\ <S |c| and 
2) \d\ ^ |c|. In order to study these cases, we assume J = or c = respectively. 
The results of our fits for the P' > {P' < 0) cases are shown in Fig. 3a (Fig.3b). 
Solid (dashed) lines correspond to J = (c = 0). 

Clarification of how the curves in Fig. 3 were obtained is in order. The ap- 
proximation leading to Eq.(22) was based on the assumption that FSIs may be 
treated perturbatively. Consequently, the FSI parameters d, c cannot be too large. 
Consider for example the dT correction to the penguin SD amplitude P. Since 
the ratio of iPj/lTl is expected to be around 0.3 (in our fits without FSIs we have 
= 0.73/2.58 = 0.28), the admissible value of \d\ should be smaller than that number. 
Consequently, in the analysis leading to Fig. 3 wc limited the region of parameter 
values to |Re d\ < 0.25, |lm d\ < 0.25. In order to give a feeling for the expected 
scale of FSI parameters, let us recall that the contribution to arising from quasi- 
elastic non-Pomcron rescattering is fully calculable in the Regge model, and in ref. [8] 
it was estimated to be of the order of 0.04 — 0.05 . The value of \d\ of the order of 0.1 
or 0.2 could therefore represent the sum of contributions from several intermediate 
channels while being still acceptable as corresponding to a perturbative realm. 

When our restrictions on the allowed values of |Re d\, |Im d\ are relaxed, the 
global minima seen in Fig. 3 are still present with the same values of /. For the 
P' > case (Fig. 3a), the relevant curve lies only slightly below that shown. For the 
P' < case (Fig. 3b), the minimum of the dashed curve on the right (at 7 ^ 130") 
becomes deeper with the vahic of / comparable to its value at ^ ^ 60°. However, 
the corresponding value of \d\ becomes significantly larger than 0.25. The fitted 
values of |c| are of the order of 0.25 also when |c| is not restricted. In the presented 
fits no restrictions on c were imposed. 

Comparison with Fig. 2 shows that the minima of / treated as a function of 7 
are now deeper and significantly shifted when compared with the no-FSI case. 
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For the P' > case, we have: in the d = case the minimum of /(7) appears 
at 7 ~ 50° with a value of / at minimum being 12.2 and c = —0.28 (Fig.3a, sohd 
hue), while in the c = case 7 ~ 80" is singled out with / = 13.3 and lie d = 
+0.25, Im d = —0.21 (Fig.3a, dashed line). The reduced Xred — f/i-^ ~ ^)) with 
= 15 used as the number of data points, and k being the number of independent 
parameters, goes down from Xred around 25/(15 — A) ^ 2.2 for the no-FSI case to 
Xred around 1.2 — 1.4 when FSI is taken into account. 

For the P' < case, the minima of /(7) are significantly deeper: in the d = 
case there is a slight shift in 7 (from around 102*^ to around 90° ) with the value 
of 7(7) at minimum being 8.84 and c — 0.24 (Fig. 3b, solid line); in the c — case 
the shift in 7 is larger and a minimum appears at 7 = 57° with the value of /(7) 
at minimum being 7.61 (Fig. 3b, dashed line). In the latter case the fitted values of 
FSI parameters are 

Red -0.22 (25) 
Imd ^ +0.21 (26) 

In both cases the value of Xred about 0.9. The second minimum of the dashed hne 
in Fig. 3b at 7 130° corresponds to a different sign of Re d. When the restriction 
on the size of \d\ is relaxed, this minimum becomes as deep as that at 7 = 57°. 
Then, however, the value of \d\ is much larger than 0.25. Since Xred significantly 
smaller for P' < 0, we restrict further discussion to this case. 

In Fig. 4, relaxing for a moment the assumption e = 0, we show the e-dependence 
of the minimal values of / for P' < and for fixed values of 7 in the two cases 
of J = (Fig.4a) and c — (Fig.4b). The region of small e seems to be preferred 
in both cases. In this analysis, as in that leading to Fig.3, the values of d were 
restricted to |Re d\ < 0.25, |Im d\ < 0.25, while the values of c were set free. 

In the most general fit (with P' < 0), we assumed e = and simultaneously 
treated all three FSI parameters (Re d, Im d, c) as free. In Fig. 5a we show the 
contour plot of the minimum of / treated as a function of complex d. The fitted 
values of c are not shown but in the region around Re d — —0.22, Im J = +0.21 
(point X) they turn out to be close to 0. Thus, allowing c to be free does not lead 
far away from the minimum found before for the c — case.The corresponding Xred 
is around 1.0. The fitted values of |c| turn out to be smaller than 0.25 for all of 
d in Fig.5a with the exception of a thin slice on the right (for Re J > 0.20 and 
lmd< 0.05). 

In order to show what happens for other negative as well as for positive Re d, 
below we present also fits performed at two additional points (pi) and (p2): 
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and 



point pi : He d = —0.10 
Imd = +0.15 

point p2: Rc(J=+0.15 
Imd^ +0.15 



(27) 



(28) 



The B PP branching ratios corresponding to the four cases (d = 0, c = 0, point 
(pi), point {p2)) are gathered in Table 3 together with other fit details. 

As can be seen from Table 3, the quality of the description of branching ratios at 
points (pi), (p2) is essentially the same as that at minimum (point X, c ^ 0). Table 
3 shows also that the dominant contribution to / comes from the 2(T discrepancy 
between the experimental and the fitted — > n^K^ branching ratios (a similar 
problem with this decay channel can be observed in other papers, see eg. [33] ). In 
a recent paper [39] , the question of a potential discrepancy in the sum rule relating 
the branching ratios in B~^,B'^ — > Ktt decays was discussed and it was suggested 
that the experiment hints at a slight enhancement of both modes involving 7r°. In 
our fits (as in [33]), however, the measured branching ratio of B^ n^K'^ is well 
described. 

Fig. 5b gives the contour plot of the corresponding fitted values of 7. In the 
region around points X and pi the fitted values of 7 seem to be in agreement with 
the conservative SM expectation of 50° < 75M < 80°, so this part of the complex d 
plane may be called the " SM" region. 



5.3 CP asymmetries 

With the values of the FSI (and other) parameters fixed, one can attempt the 
calculation of CP-violating observables. The CP-violating asymmetries in S — > Kn 
decays defined as 

{B — B^,B^, B — B^,B~) were calculated for all four cases under discussion. 
The relevant predictions are given in Table 4 together with the experimental data 
([31, 40, 41, 42, 43, 44, 45]) as averaged in [46]. The "SM" region of small c and 
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negative Re d (represented by points X {c — 0) and pi (c = —0.11)) seems to 
describe the experimental B — > Ktt CP asymmetries somewhat better than the 
d = case or the region of positive Re d (ie. point p2) do: our FSI approach 
prefers negative 5° tt~K^ asymmetry, in agreement with the experiment and in 
disagreement with the predictions of ref.[l]. Although the B Kit asymmetries 
are experimentally small, they might provide important model tests (see eg. [11]). 

In view of the recent BaBar measurement [47] favouring a large negative asym- 
metry in B^ TT'^r] decays, we have computed the asymmetries in all B~^ 7r~^7]{ri') 
and B~^ K~^r][r]') decays. The results are given in Table 5 with the data [47, 48, 49] 
averaged as in [46]. Contrary to the BaBar result, our B~^ — > 7i~^r} asymmetry is 
small and positive for 7 in the " SM" region. On the other hand, our K~^r] asymme- 
try (fairly large when compared with other asymmetries) seems to agree with the 
data. Problems with the simultaneous description of 7r~^r] and K~^r) asymmetries 
have been noted in [46] as well. 



We have also calculated parameters relevant for the time-dependent rates in 
^ JT+x- [50], ie.: 

i -+- IAttttI 

and 

1 _ |A^^p 

^'^'^ ^ 1 I |\ 12' (31) 

where 

Our predictions are given in Table 6. Although the experimental results from Belle 
and BaBar [51, 52] still exhibit the well-known discrepancies [53, 54], the "SM" 
region of small (negative) c and negative Re d (with the value of 7 close to the SM 
expectations) seems favoured again. 

For the time-dependent rates in B^(t) — > i]'Ks, the effect of final-state interac- 
tions is negligible. Indeed, the relevant amplitudes are dominated by the P' and 
S' amplitudes (in particular, the FSI correction is dominated by terms proportional 
to P', see Table 2). Thus, all important terms have the same weak phase. Conse- 
quently, one obtains S^j/Rs ~ sin 2/3, Cj^'Ks ~ 0, in agreement with the experimental 
averages (from [46]) of S^>Ks = +0-33 ± 0.25, = -0.18 ± 0.16 . 

The B~^ vr+TT^ asymmetry is predicted to be zero (cf. Tables 1 and 2), in 
agreement with its experimental value of — 0.07± 0.14 (average from [46]). 

Although apart from the discrepancy in sign with the most recent BaBar ir'^r] 
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result there seems to be a hint of agreement with other asymmetries, one has to 
remember that these (and other) predictions for asymmetries may be affected by 
the inclusion of the charming penguin contribution [23, 34]. 

6 Conclusions 

In this paper we have analysed the contributions from both clastic and inelastic 
SU(3)-breaking final-state interactions in B decays to two light pseudoscalar mesons 
{B PP). 

We have found that the inclusion of an experimentally determined pattern of 
SU(3) breaking in Pomeron- induced rescattering enhances the value of the S' / P' ra- 
tio when extracting it from the fit to the B — > PP branching ratios. However, taking 
this rescattering into account does not lead to any significant change in the overall 
fit. Since at the energy of s = m\ the inclusion of non- Pomeron clastic rescattering 
may lead to small corrections only, analyses incorporating full clastic FSIs can lead 
neither to a significant improvement in the quality of data description, nor to the 
extracted value of 7 being substantially shifted towards the SM expectation. 

We have pointed out that a small value of the B — > K^K~ branching ratio does 
not imply neghgible inelastic rescattering effects in other B — > PP decays. This 
conclusion follows from the fact that rescattering in the B — > K^K~ decay is inde- 
pendent of two of the three parameters describing the totality of inelastic FSIs: one 
related to the u-loop long-distance penguin (in a resonance channel), and the other 
one describing quark rearrangement (in an exotic channel). As for B^ K~^K^, 
with U-spin symmetry probably broken by final-state interactions, this decay was 
argued to be less helpful in the determination of the size of rescattering effects than 
originally suspected. Its importance in the determination of the size of rescattering 
effects (ie. the size of the w-loop long-distance penguin) would then lie not just in 
its relation to 5+ — > n'^K^, but, more properly, in its relation to all other B — > PP 
decay channels. 

Finally, after neglecting the relative strong phases of short-distance amplitudes, 
we have carried out fits to the available B — > PP branching ratios with all elastic and 
inelastic SU(3)-breaking rescattering effects taken into account. The only neglected 
but potentially important corrections were those due to the intermediate states 
composed of charmed mesons. Our fits show the importance of rescattering effects 
and weakly hint at the value of 7 compatible with SM expectations. However, other 
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values of 7 are also possible. Narrowing the range of admitted values of 7 will require 
taking into account the experimental data on asymmetries in addition to those on 
branching ratios. In this paper we used the values of rescattering parameters as 
determined from the fits to the branching ratios, and predicted several CP-violating 
observables (CP asymmetries in B Kn decays, Sj^t^ and Cj^j^ for the B^(t) tt+tt^ 
time-dependent decay rates etc.). Again, weak agreement with the data (with the 
notable exception of the tt^t} asymmetry) was found for 7 close to the SM 

expectations. 

This work was supported in part by the Polish State Committee for Scientific 
Research (KBN) as a research project in years 2003-2006 (grant 2 P03B 046 25). 
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FIGURE CAPTIONS 



Fig. 1. Types of rescattering diagrams: (u) uncrossed, (c) crossed. 

Fig. 2. Dependence of minimized function / (Eq.(24)) on 7: thin lines - P' > 0, 
thick hues - P' < 0; sohd Uncs - no FSI/SU(3) symmetric Pomeron-induced FSI, 
dashed hues - SU(3)-breaking Pomeron-induced FSI. 

Fig. 3. Dependence of minimized function / (Eq.(24)) on 7 for full FSI: (a) 
P' > 0, (b) P' < 0; sohd hues - d — 0, unrestricted |c|; dashed hues - c — 0, 
|Re d\ < 0.25, |Im d\ < 0.25. 

Fig. 4. Dependence of minimized function / on e for full FSI and different values 
of 7: (a) J = 0, unrestricted |c|; (b) c = 0, |Re d\ < 0.25, |Im d\ < 0.25. 

Fig. 5. (a) Contour plot of minimized function / in complex d plane. Positions 
of the minimum (X) and of the selected points pi, p2 are indicated, (b) Contour 
plot of fitted values of 7 in complex d plane. 
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Table 3: Branching ratios of B decays (in units of 10 ^) 





decay 


expt 






P' < 












no FSI 


Pomeron 


d = 


c = 


pi 


p2 


B+ 


— >■ TT+vr'^ 


5.8 ± 1.0 


4.85 


4.79 


5.23 


5.38 


5.54 


5.86 






0.0 ±2.0 


0.57 


0.51 


0.54 


1.09 


1.02 


0.87 






2.9± 1.1 


2.13 


2.13 


3.47 


2.90 


2.60 


2.50 




TT'^T]' 


0.0 ±7.0 


1.06 


1.03 


1.69 


1.39 


1.25 


1.22 




— > TT'^TT" 


4.7 ±0.5 


4.93 


4.93 


5.19 


4.77 


4.79 


4.62 






1.9 ±0.7 


0.55 


0.56 


1.98 


1.85 


0.82 


1.31 




K+K- 


0.0 ±0.6 
























0.0 ±4.1 


0.53 


0.48 


0.50 


1.02 


0.95 


0.87 


B+ 




18.1 ± 1.7 


18.28 


18.51 


19.70 


19.15 


18.98 


20.53 






12.7 ± 1.2 


12.96 


12.87 


12.47 


12.15 


12.34 


12.76 




r]K+ 


4.1 ± 1.1 


2.45 


3.05 


3.64 


4.18 


4.07 


4.24 




7]'K+ 


75.0 ±7.0 


72.85 


72.09 


69.31 


69.07 


69.53 


69.60 




7r-K+ 


18.5 ± 1.0 


18.90 


18.90 


17.57 


18.89 


18.99 


18.10 






10.2 ± 1.5 


6.38 


6.53 


6.79 


7.16 


7.04 


7.37 






0.0 ±9.3 


1.83 


2.43 


4.28 


2.50 


2.29 


5.36 






56.0 ±9.0 


67.07 


66.62 


65.68 


66.51 


65.37 


65.06 




fmin 




16.05 


14.25 


8.84 


7.61 


9.70 


8.86 




\T\ 




2.58 


2.56 


2.41 


2.71 


2.69 


2.66 




P' 




-4.14 


-4.24 


-4.34 


-6.17 


-5.98 


-5.53 




S' 




-1.77 


-2.27 


-2.09 


-1.53 


-1.41 


-1.52 




Ifit 




103° 


101° 


89° 


57° 


78° 


99° 




c 








±0.24 





-0.11 


±0.18 




Re d 











-0.22 


-0.10 


±0.15 




Im d 











±0.21 


±0.15 


±0.15 



Table 4: Asymmetries in — > Kir decays 



decay expt 


P' < 

J=0 c = pi p2 


B+ n+K^ -0.032 ± 0.066 
B+ Tr°K+ ±0.035 ± 0.071 
5° ^ 7r-K+ -0.088 ± 0.040 
S° ^ 7r°i^° 0.03 ± 0.37 


±0.09 ±0.05 -0.07 
-0.04 -0.10 -0.03 ±0.03 
±0.03 -0.10 -0.07 ±0.08 
±0.07 ±0.13 ±0.04 -0.05 
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Table 5: Asymmetries in tt"*" 77(77') and K~^r){r)') decays 



decay expt 


P' < 

J=0 c = pi p2 


B+^TT+T] -0.51 ±0.19 
B+ n+ri' 

B+ K+r] -0.32 ± 0.20 
B+ K+ri' -0.002 ± 0.040 


+0.10 +0.06 -0.09 
+0.10 +0.06 -0.10 
+0.23 -0.39 -0.49 +0.32 
-0.01 +0.01 +0.01 -0.01 



Table 6: CP-violating parameters in time-dependent rates for B — > tt+tt 



parameter experiment 
Belle 
BaBar 


P' <0 

J=0 c = pi p2 


S^^ -1.23 + 0.411°:°? 

-0.40 + 0.22 + 0.33 
C^^ -0.77 + 0.27 + 0.08 

-0.19 + 0.19 + 0.05 


-0.12 -0.78 -0.23 +0.49 
-0.05 -0.21 -0.08 +0.11 
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Figure 5: 
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